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Abstract 

We present a dynamical model of lipoprotein metabolism derived by combining a 
cascading process in the blood stream and cellular level regulatory dynamics. We 
analyse the existence and stability of equilibria and show that this low-dimensional, 
nonlinear model exhibits bistability between a low and a high cholesterol state. A 
sensitivity analysis indicates that the intracellular concentration of cholesterol is 
robust to parametric variations while the plasma cholesterol can vary widely. We 
show how the dynamical response to time-dependent inputs can be used to diagnose 
the state of the system. We also establish the connection between parameters in the 
system and medical and genetic conditions. 

Key words: lipoprotein metabolism, dynamical systems, nonlinear models, 
metabolic control mechanisms, statins, familial hypercholesterolemia 



1 Introduction 



We model the dynamics of lipoprotein metabolism associated with the trans- 
port of lipids in the blood. The malfunction of the delivery of lipids and 
cholesterol from the liver to the cells has important ramifications. There is 
strong medical evidence linking high plasma concentrations of low-density 
lipoproteins (LDL) to the d evelopment of atherosclerosis, the deadliest disease 
in industrialised countries ([Glass and Witztuml . 12001c iRodriguez et all Il999c 



LibbvllMav 2 002). However, LDL is just one species in a metabolic cascade by 



which different lipoproteins (Table 1) are synthesised, degraded and absorbed. 
In this complex network, it is important to investigate how metabolic control 
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is exercised and if dynamical features, other than high LDL concentrations, 
could be useful markers for the detection of disease states. 



Previous modelling efforts of lipoprotein metabolism vary in scope. Some mod- 
els have concentrated on detailed aspects of the syst em, such as the fluid dy- 
namics of lipid accumulation on the arteri al walls ([Hazel and Pedlev . Il998h 
or the chemical kinetics of LDL oxidation ( Cobbold et all l2002j ). but do not 
describe the processes at the cellular or physiological level. Other approaches 
have model led the l i poprotein networ k as a whole with linear co mpartmen- 
tal models (|,Tacaue A flQsi IPont et all Il998l : IParhofer et all Il99lh . Although 
these models can be fitted to match experimental data, the compartments lack 
rigorous physiological meaning, thus making it difficult to describe the under- 
lying biochemical processes. For instance, compartments added to represent 
subclasses of lipoproteins can always be reduced to one com partment because 
of the linearity of the models ((Eisenfled and Grundvl Il984l ). 



Our model of lipoprotein metabolism is a relatively low-dimensional system 
of nonlinear differential equations which are strongly linked to the underlying 
physiological processes. The system can undergo a transition between low 
and high LDL steady states, a feature which is robust to randomness and 
uncertainty in the parameters. Because the equations are directly related to 
physiological processes, the model can be used to study the effect of genetic 
or behavioral conditions and could serve as an aid to test hypotheses for 
diagnosis and intervention. In what follows, we derive our model from the main 
physiological processes; we then analyse the system and carry out a sensitivity 
analysis; finally, we show how the results can be related to experiments and 
medical conditions. 



2 Modelling Lipoprotein Metabolism 



Lipoproteins constitute the primary means of transport of lipids from the liver 
to the cells via the blood stream. Lipoprotein dynamics is intimately connected 
to lipid and cholesterol metabolic networks. Figure 1 provides a schematic 
overview of the origin, transport and fate of the different lipoproteins in the 
human body at two levels: the macroscopic and the cellular. 

Lipoproteins are globular aggregates of lipids and proteins: cholesterol ester 
and triacylglycerol form the core; amphiphilic phospholipids surround them 
forming the hull; and different apolipoproteins (apo A, apo B, apo C, and apo 
E), which function as trigger-molecules for specific reactions, are embedded 
in the surface. Lipoproteins are classified in five standard groups (Table 1), 
based on density and the apolipoproteins attached to them. In increasing or- 
der of density and decreasing order of size (the fewer lipids relative to its size, 
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the denser the lipoprotein is): chylomicrons (Chyl), very low-density lipopro- 
teins (VLDL), intermediate-density lipoproteins (IDL), low-density lipopro- 
teins (LDL), and high-density lipoprotein (HDL). 





Chyl 


VLDL 


IDL 


LDL 


HDL 


Density (g/ml) 


< 0.95 


0.95-1.006 


1.006-1.019 


1.019-1.063 


1.063-1.21 


Diameter (nm) 


80-100 


30-80 


25-30 


20-25 


8-13 


Cholesterol (%) 


3-5 


12-21 


27-46 


40-50 


15-25 


Apolipoproteins 


A, C, E, B-48 


C, E, B-100 


C, E, B-100 


B-100 


A, C, E 



Table 1 

The properties and composition of the five standard lipoprotein 
classes (|AdielsL 12002) : chylomicrons (Chyl), very low-density lipoproteins (VLDL), 
intermediate-density lipoproteins (IDL), low-density lipoproteins (LDL), and high- 
density lipoproteins (HDL). 



Chylomicrons contain a large proportion of fatty acids and are found in the 
blood stream mainly after digestion of a meal. They are synthesised by the 
intestines and constitute the direct means by which dietary fat is delivered to 
heart, muscle, and adipose tissue. Chylomicrons provide the liver with lipids 
and proteins but are not regulated by the liver. Therefore, we do not consider 
the concentration of chylomicrons as an explicit variable in our model and 
include their effect parametrically as an input of lipids to the system. 

The liver is responsible for the secretion of all other lipoproteins, most sig- 
nificantly VLDL. (Although small amounts of IDL and LDL are produced by 
the liver, we neglect this contribution.) VLDL is secreted by the liver into 
the blood at an almost constant rate (slightly higher during the night) as a 
means of transport of synthesised fat and cholesterol to peripheral tissues. 
While in the blood stream, the apo C protein on the surface of VLDL can 
bind to the extracellular enzyme lipoprotein lipase (LPL), which is present 
on the capillary wall. The effect of LPL is to hydrolyse triacylglycerol, thus 
reducing the size of VLDL and increasing the percentage of denser molecules 
(i.e., it degrades VLDL to IDL). After hydrolysis, the fatty acids are taken up 
by nearby cells or by serum albumin for transport to more peripheral cells. A 
similar process leads from IDL to LDL. This is the lipoprotein cascade that 
transforms VLDL into IDL into LDL. (This process also involves losing apo 
C and apo A to HDL, which is independently synthesised in the liver and is 
not a part of the lipoprotein cascade.) 

About 50% of the IDL and 75% of LDL in the blood are absorbed by hepatic 
and peripheral cells via LDL receptors (LR) on the cell membrane. These re- 
ceptors are synthesised by the cell and recognise apo B and apo E with high 
affinity. In liver cells, the absorbed LDL is reused for lipoprotein synthesis and 
excess cholesterol is secreted into bile. In non-hepatic cells, the absorbed LDL 
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Fig. 1. Overview of origin, transport and fate of lipoproteins. A schematic 
view of the processes involved in the metabolism of lipoproteins at both the macro- 
scopic and cellular levels. The liver secretes VLDL at a rate uy. The turnover rates 
of VLDL and IDL in the blood stream are ky and hi, respectively. At the cellular 
level, d,L and dj denote rate constants of the receptor-mediated uptake of lipopro- 
teins: IDL and LDL bind to membrane receptors (LR), which cluster in coated pits 
in the cell membrane. Note that it is also possible for LDL (but not for IDL) to be 
absorbed directly through a non-receptor mediated path with rate constant d. After 
endocytic vesicles are formed and internalised, the lipoproteins are hydrolysed in the 
cell releasing lipids into the cytoplasm. A significant proportion of these lipids (xi 
in IDL and xl m LDL) is cholesterol which contributes to the level of intracellular 
cholesterol (IC). The intracellular cholesterol is used for cell function or eliminated, 
mainly through the action of HDL, at a rate die- Most of the LDL receptors (but 
not all) are reincorporated into the membrane. In addition, the nuclear synthesis 
of the receptors LR is regulated through negative feedback (h) by the intracellular 
cholesterol IC. 
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supplies the cholesterol content esse ntial for cell function ((White and Baxter 
1984 : Chopra and Thurnham . 1999f ). Although almost every cell can synthe- 
sise cholesterol to some extent, we assume that all cholesterol has to be deliv- 



ered to_y^_ceUs_j2a_JjDL_or IDL absorption (Mathe ws et all l2000t iCooper 
2000UConverse and Skinnerl . Il99i llDietschv et all Il978h . 



As mentioned above, there is another important component in the lipoprotein 
network: HDL (the "good cholesterol"), which is also secreted by the liver. 
HDL is responsible for reverse cholesterol transport, the transport of excess 
cholesterol from cells and other lipoproteins back to the liver. The effect of 
HDL is included in our model through a rate constant that parameterises the 
elimination of intracellular cholesterol, which is significantly affected by the 
HDL pathway. 

The above processes can be modelled with the following system of differen- 
tial equations, where [•] denotes concentration and (pLR is the dimensionless 
fraction of total possible LDL receptors: 



d[VLDL) 

dt 
d[IDL) 

Jt 
d[LDL] 



-k v [VLDL}+u v 

k v [VLDL] - ki [IDL] - d I [IDL](f) IjR 



dt 



ki [IDL] - d L [LDL]<p LR - d[LDL] 



>LR 



dt 
d[IC] 



-b(d I [IDL]+d L [LDL])4 >LR + c 



1 - 



>LR 



[IC] 



(1) 

(2) 
(3) 
(4) 



dt 



[Xidi[IDL] + X Ld L [LDL]) $ LR + XL d[LDL] - d IC [IC]. (5) 



The linear terms ky[VLDL] and ki[IDL] in (l)-(3) correspond to the cascade 
degradation of VLDL to IDL to LDL assuming saturated kinetics (i.e., excess 
of LPL enzyme) for the reactions. An implicit assumption is that the blood 
flow does not play a significant role in the absorption process. This is based 
on our analysis of the relevant time scales and also on detailed calcu lations of 
mass transfer in tissue under physiological conditions (Zhang, 20051 ). 



Equation (1) contains the only input to the system, uy, which represents the 
overall secretion rate of VLDL by t he liver. High value s of uy are associated 



with a high dietary intake of fats ( Spadv et al. . 1993| ). On the other hand 



medication with statins would reduce uy by lowering the synthesis of choles- 
terol and VLDL in the liver. 



The nonlinear terms dj[IDL]<pL R and dL[LDL](f)L R represent the endocytosis 
of IDL and LDL via a process involving LDL receptors (LR) and leading to 
the release of lipids and cholesterol into the cytoplasm. This c an be modelled 
in its simplest, classical form through a second order reaction ((Bvwater et al. . 
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20021 ). In addition to this receptor- mediated uptake, mammalian cells can take 
up LDL (but not IDL) by non-specific endocytosis at a ra te d\LDL], linearh 



up lul [but not lul) Dy non-specmc enaocytosis at a ra te a\LUL\, linearly 
proportional to its plasma (or extracellula r ) concentration dHobbs et a . , 1990; 



Willnow et all Il999l : iQlass and Witztuml . l200lt iQoldstein and Brownl . I'l997) . 



Endocytosis is greatly reduced under the genetic disease familial hypercholes- 
terolemia (Section 4.1.3). 

The last two equations (4)— (5) describe cholesterol uptake and regulation at 
the cellular level. The variables are: 4>lr, the fraction of the total possible LDL 
receptors that mediate IDL and LDL endocytoses; and [IC], the concentration 
of cellular cholesterol, which closes the loop by acting as a cellular control for 
the expression and synthesis of the LDL receptors, LR. 



Parameter 


Source 


Units 


Nominal Value 


Range 


ky 


(Packard et al.. 2000) 


h- 1 


0.3 


0.15-0.6 


kj 


(Packard et al.. 2000) 


h- 1 


0.1 


0.025-0.1 


dj 




h- 1 


1.4 


0.5-2 


d L 




h- 1 


0.005 


0.005-0.02 


d 


(Dietschv et al.. 1993) 


h- 1 


0.0025 


0.0025-0.0075 


b 




ig- 1 


0.1 


0-2 


c 


( Goldstein and Brownl. Il977l) 


g{lh)- 1 


0.05 


0-1 


XI 


(Aciiels, 20j]2) 


0.35 


0.25-0.45 


XL 


(Adiels. 2002) 




0.45 


0.4-0.5 


Uy 


(White and Baxter. 1984) 


g{lh)- 1 


0.3 


Variable 


die 


(White and Baxter. 1984) 


h- 1 


0.45 


Variable 



Table 2 

Parameters of the model (l)-(5). Nominal values and ranges as found in the 
literature. The nominal values for di and di are obtained by assuming that all other 
parameters are at nominal values and the equilibrium point of t he linearised system 
in th e lower branch of Eq. (17) is that of normal individuals ([White and Baxter . 
1984): [VLDL] = lgl' 1 , [IDL] = 0.2 gl' 1 , [LDL] = 2gl~ 1 . The nominal values 
for uy and die are crude estimates from the cited reference. Because of their high 
dependence on diet, medication and genetic factors, we consider uy and die as the 
control parameters. 



Eq. (4) reflects the dynamical processes involving the LDL surface receptors, 
which take part in the endocytotic reactions. Experiments in normal individ- 
uals indicate that the maxim um number of receptors on th e surface of the cell 
is 15,000 — 70,000 per cell ([Goldstein and Brownl . Il977l ). The variable 4>lr 
represents the fraction of the total possible surface receptors that are present 
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on the surface of the cell. Therefore, 4>lr is a surface occupancy fraction; hence 
it is dimensionless and bounded (0 < 4>lr < 1). Biological evidence indicates 
that most (but not all) of the receptors used in the endocytosis are reinte- 
grated into the cell membrane. The first term in (4), weighted by a factor 
b, represents the loss of LDL receptors that are not recycled. The last term 
in (4) describes the cholesterol-regulated synthesis of LDL receptors and the 
process of their attachment to the surface of the cell. Further biological evi- 
dence indicates that cells regulate the synthesis of LDL receptors to prevent 
the over-accumulation of intracellular cholesterol. When grown in the presence 
of different concentrations of LDL, cells adjust the number of LDL receptors to 
take up o nly enough LDL to sat i sfy th e cholesterol requirement for membrane 



turnover ([Goldstein and Brownl . 119771 ): when grown in total absence of LDL, 



the maximum number of possible surface LDL receptors is reached after 2-3 
days. We have modelled this inhibitory mechanism with a simple reciprocal 
feedback: the synthesis of cytosolic LDL receptors is inversely proportional 
to the intracellular cholesterol [IC]. Furthermore, the attachment of the cy- 
tosolic LDL receptors to the cell surface can be described in its simplest form 
with a rate that is proportional to the fraction of unoccupied receptor sites 
(1 — (Plr)- The combined rate at which (ftm increases is then given by the 
product c(l — (pLji)/[IC], where c modulates the weight of the combined pro- 
cess of regulated synthesis and attachment. However, we emphasise that the 
qualitative results of the model do not depend on the particular form of the 
feedback chosen, as we discuss at the end of Section 3.2. For a more detailed 
explanation of this term see Appendix A. 

Finally, Eq. (5) establishes the balance of flows of the cellular cholesterol [IC] . 
The parameters < xi < Xl < 1 are the proportions of cholesterol in IDL and 
LDL, respectively, which can be measured experimentally. The last term in the 
equation djc[IC] represents the outflow rate of cholesterol from the cell. This 
is largely related to the HDL-mediated reverse cholesterol transport, but also 
accounts for the cholesterol used by cells for their metabolism or eliminated 
through bile secretion in the liver. Therefore, our model parameterises the 
effect of HDL through the constant die- As in the case of uy, die will be a 
highly variable control parameter which can be affected by diet and statins. 

The model has eleven parameters with nominal values and ranges shown in 
Table 2: 



Model parameters: k v , ki, di, di, d, b, c, xi-, an d Xl- These will differ 
between individuals but are assumed to remain virtually unchanged for an 
individual except for slow changes (e.g., with age). We have found values 
and estimates in the literature for all of these parameters, except for b for 
which we assume b = 0.1. Sensitivity analyses will be performed to check 
that the qualitative results are robust to variations in these parameters. 

Control parameters: uy and djc- The secretion rate of VLDL from the 
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liver (uy) and the depletion rate of intracellular cholesterol {die) are highly 
variable parameters and affected by genetic factors, diet and medication. 
Hence, we will consider these as the control parameters for our analysis. 



3 Model Analysis: Fixed Points and Sensitivity Analysis 



3.1 The System at Equilibrium 



We have carried out a numerical analysis of the bifurcations and fixed points 
of the dynamical system (l)-(5). Before presenting the numerics, we briefly 
note some analytical features of the model. 

The system is effectively four- dimensional since [VLDL] converges exponen- 
tially to its equilibrium value [VLDL] eq = uy/ky. Moreover, we can restrict 
our analysis to the non-negative orthant because the system is essentially 
non-ne gative, i.e., solutions remain non-negative for non-negative initial con- 
ditions ( Haddad et al. . 200lh . The equilibrium point can also be bounded from 



above using the nullclines of (l)-(5) to give: 



klUV <[LDL] eq <^ (7) 



(fcj + di){d + d L ) iq ~ d 

0< <f>LR eq <1 (8) 
Xidi + XhkiUy XhUy , . 

k I + dj dTc - [ICU - ^ (9) 

The equilibrium approaches the bounds (6)-(9) asymptotically as uy — > and 
uy — > oo, as shown in Figure 2. The simple linear bounds (6)-(9) represent well 
the asymptotic behaviour although tighter nonlinear bounds (not shown) can 
be obtained. In conclusion, there is a compact, connected, positively invariant 
bounding region for the system. This implies that any linearly stable solution 
in this region is also locally asymptotically stable. 

Consider now the fixed points of system (l)-(5) which are the solutions to the 
following algebraic equation: 

^LR eq diccd 2 jd L + (p 3 L R eq (d Ic cd I (2k I d L + di{d - d L )) + Uybd]d L xi) 
+ (j> 2 LReq (dicc(2kidi(d - d L ) + kjd L - djd) + u v bd I {d L ki(xi + Xl) + ddiXi)) 
+ (f>LR eq {di C cki{ki(d - d L ) - 2did) + u v bk I (dd I XL + cfcfc/Xz)) - d Ic ckjd = 

(10) 
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together with 



\VLDL\ 



eg 



Uy 



\IC) 



ki + di4> LRi 
( X idj[IDL] eq + X Ld L [LDL) eq 



[LDL] 



kj[IDL] 



(:•(/ 



(:'</ 



+ X Ld[LDL] eq 



eq 



die 



(11) 



Descartes' rule of signs implies that the number of positive equilibria is 1 
or 3. If we assume that the system has neither (quasi)periodic solutions nor 
strange attractors (which never appear in our numerics), then if there is only 
one fixed point, it is globally asymptotically stable; and if there are three, two 
are locally asymptotically stable and one is unstable and lies "between" the 
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Fig. 2. Plot of equilibrium solutions versus uy. We show stable solutions 
(solid lines), unstable solutions (dotted lines), and bounds (6)-(9) (dashed lines) 
as a function of uy with parameters: { ky = 0.3/i -1 , kj = 0.025/i _1 , dj = 2h~ 1 , 
d L = 0.01/T 1 , d = 0.0075/1- 1 , xi = 0.1, XL = 0.6, b = 0.1/cT 1 , c = 0.05c/(Z/i)- 1 ; 
die = lh }■ For this set of parameters the system shows bistability. The bifur- 
cations, obtained through (14), occur at: u^bif ~ 1-637 and uy,bif ~ 1-811. The 
approximation (16) to the latter is 1.759. 

We numerically explore the fixed point equation (10) and record how the num- 
ber and stability of positive fixed points change as the control parameters are 
varied. Figure 2 shows a numerical bifurcation diagram as a function of uy, 
for a given set of parameters. The system has two distinct (sometimes coex- 
isting) stable solutions: a low and a high cholesterol branch. For small (high) 
uy, only the low (high) branch exists. For intermediate values of Uy, the two 
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stable branches (high and low) and an unstable solution coexist. The existence 
of this coexistence region leads to the possibility of hysteresis and a delayed 
return to the normal state when uy is swept. The transitions that delimit 
the coexistence region are two saddle-node bifurcations. A similar bifurcation 
diagram (Figure 3), with equivalent saddle-node bifurcations, is obtained by 
sweeping the other control parameter die, although in this case the transition 
to a higher cholesterol state is produced by decreasing die- These bifurcations 
could be of physiological relevance because the parameters uy and d IC can be 
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d, c (h- 1 ) 



d, c (h- 1 ) 



1.5 



1.5 



Fig. 3. Equilibrium solutions versus die- We show stable solutions (solid lines), 
unstable solutions (dotted lines), and bounds (dashed lines) as a function of die- 
All parameters as in Fig. 2 and uy = Igilh)" 1 . The bifurcations for this set of 
parameters, obtained through (14), occur at: djc.ba ~ 0.373 and d/cbif ~ 0.305. 
The approximation (16) to the latter is 0.323. 

The nature of the bifurcation diagram depends on the non-control parame- 
ters of the model. In fact, the coexistence region does not exist in certain 
regions of parameter space. In that case, the transition from the low to the 
high cholesterol solution is smooth, without undergoing a bifurcation (i.e., no 
sudden jumps and no hysteresis). We make this more precise by analysing the 
equation for the bifurcation points. Rewrite the fixed point equation (10) as: 



(y - \){y + g){y + ff + Z y (y + hf)(y + f + g) = 



(12) 



where y = 4>LR eq an d the following ratios of parameters have been defined: 



e = b/c, f = h/di, g = d/d L , h = xl/xi, Z = u v exi/d IC - 



(13) 
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The bifurcation points are the solutions of an augmented system of equations, 
formed by the fixed point equation and the condition that the Jacobian of (1)- 
(5) vanishes. Using Groebner bases (or considerations about monotonicity and 
curvature of (12)) we obtain a quintic polynomial equation for the bifurcation 
points y hi f- 



y ^ + (2g + 2fh + f)yt ii +(g 2 + 3f 2 h + f-fh-gf + 4fhg)yl i{ 
+ ^gf + f 2 -f 2 g + fh-2fhg + 2f 2 hg-g 2 f-f 2 h + 2fhg 2 )y 2 hi{ 
+ (2f 2 g-fhg 2 + 2g 2 f + f 2 hg) y m + f h g + f 2 h g 2 = 0. (14) 



Based on Descartes' rule, we conclude that there is no coexistence region (and 
no bifurcations) if: 



gi + 3f 2 h + f-fh-gf + 4fhg>0 
4gf + f 2 -f 2 g + f 3 h-2fhg + 2f 2 hg-g 2 f-f 2 h + 2fhg 2 >0 

2f 2 g-fhg 2 + 2g 2 f + f 2 hg>0. 



These inequalities can be tested for each set of parameters to establish a 
sufficient condition for the saddle-node bifurcations to exist. 

Further insight can be gathered through explicit approximate formulas for the 
bifurcations. Consider the limit case with al = 0, i.e., when there is no direct 
uptake of cholesterol by the cell and only receptor-mediated intake is possible. 
In this case, the fixed point equation (12) simplifies to: 



Eq. (16) marks the point above which the low cholesterol branch ceases to 
exist and serves as a good approximation for the bifurcation of the full model if 
d is small. The approximate formulas (15)— (16) explain the overall dependence 
on uy and die of the bifurcation points shown in Figs. 2 and 3. 

3.2 Sensitivity Analysis of the Equilibrium Solution 

In order to ascertain how variations in parameter values (e.g., among individ- 
uals or due to age) affect the results, we perform a sensitivity analysis of the 



y(y + f) (y 2 + (f-i + z) y + f(Zh-i))=o, 



leading to the following bifurcation points for this limit case: 




(15) 



f(Zh-l) 
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Parameter 


Variable 


Sensitivity over lower branch (uy < "Uy.max) 


Median 


Range 


ki G [0.025,0.1] 

(Wy.max = 1-1) 


[IDL] 
[LDL] 

4>LR 

[IC] 


0.027 
1.18 

0.013 
0.17 


[0.0038, 0.036] 
[1.16,1.26] 
[3.1 x 10- 4 , 0.074] 
[0.15,0.21] 


di G [0.5,2] 

(%,max = 1-4) 


[IDL] 
[LDL] 

4>LR 

[IC] 


1.83 
1.85 
0.023 
0.18 


[1.79,2.23] 
[1.79,2.48] 
[3.1 x 1(T 4 ,0.22] 
[0.16,0.34] 


d L G [0.005,0.02] 

("V.max = 1-3) 


[IDL] 
[LDL] 

4>LR 

[IC] 


5.3 x 10~ 4 
0.82 

5.4 x 10~ 4 
3.4 x 10~ 5 


[8.5 x 10" 6 , 0.0038] 

[0.67, 0.87] 
[8.6 x 10~ 6 , 0.0039] 
[5.0 x 10~ 7 ,3.4 x 10~ 4 ] 


d G [0.0025,0.0075] 

(«y,max = I- 4 ) 


[IDL] 
[LDL] 

4>LR 

[IC] 


4.6 x 10~ 4 
0.38 

4.7 x 10~ 4 
3.0 x 10~ 5 


[5.9 x 10 -6 , 0.0045] 

[0.34,0.50] 
[6.0 x 10- 6 , 0.0046] 
[3.5 x 10~ 7 ,4.3 x 10~ 4 ] 


e = 6/cG [0.25,2.5] 
(My, m ax = !-4) 


[IDL] 
[LDL] 

4>LR 

[IC] 


0.15 
0.23 
0.15 
0.0093 


[0.0024,0.80] 
[0.0037, 1.23] 
[0.0024,0.70] 
[1.4 x 10" 4 , 0.064] 


X/G [0.25,0.45] 

(wy,max = 0.9) 


[IDL] 
[LDL] 

<f>LR 

[IC] 


0.12 
0.18 
0.12 
0.56 


[0.0039,0.76] 
[0.0062, 1.06] 
[0.004,0.71] 
[0.56,0.56] 


XL € [0.4,0.5] 

(My, m ax = 1-4) 


[IDL] 
[LDL] 

4>LR 

[IC] 


0.0018 
0.0027 
0.0018 
0.013 


[2.4 x 10" 5 , 0.014] 
[3.8 x 10- 5 ,0.21] 
[2.5 x 10~ 5 , 0.015] 
[0.012,0.021] 



Table 3 



Sensitivity of the variables of the model on the lower branch to variation 
of parameters one at a time. The sensitivity is calculated as the range divided 
by the median of each variable for five different values of the parameter that appears 
in the first column. The third and fourth columns report the median sensitivity and 
the range of calculated sensitivities over the lower branch (up to the value u^ max ). 
The salient feature of the analysis is the high sensitivity of [LDL] to all parameters 
and the relative insensitivity of [IC] to almost all parameters. 
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system. Because of its physiological relevance, we focus on the variability of 
the low cholesterol branch under parameter uncertainty. 



Table 3 summarises the sensitivity of the low cholesterol branch when pa- 
rameters are modified one at a time. Consider for instance the effect of d, 
which parameterises the non-specific uptake of LDL. We have calculated nu- 
merical bifurcation diagrams similar to those in Figure 2 for values of d e 
[0.0025,0.0075]. As d is varied, the values of the lower and upper branches 
and the location of the bifurcations change. However, as expected from (6)-(9), 
only the [LDL] is affected significantly. To make this more precise, we calcu- 
late the sensitivity for each variable over the lower branch (uy < l-4<7(Z/i) -1 ) 
by obtaining the average ratio of the range and median of the equilibrium val- 
ues. Clearly, [LDL] has the highest sensitivity to d. The sensitivity to other 
parameters is summarised in Table 3: d^ virtually only affects [LDL], as in- 
dicated by (6)-(9); variations of e = b/c displace the bifurcation point (16), 
but [IC] is fairly insensitive; variations in fcj and di have a large effect on 
[LDL], but only dj affects [IDL]\ the sensitivity of [IC] to xi is moderate as 
compared to the sensitivity of [LDL] or [IDL] to kj and dj. 

We have also performed a Monte Carlo sensitivity analysis with respect to 
variations of all parameters at once by sampling and averaging over the uni- 
form hypercube of non-control parameters in Table 3. The results show that 
the median sensitivity of [LDL] to global parametric variations is an order of 
magnitude larger than that of [IC] over the low cholesterol branch. The con- 
clusion of the sensitivity analysis is that while the plasma cholesterol [LDL] 
is sensitive to most parameters, the intracellular cholesterol [IC] is robust to 
parameter perturbations. In essence, this is a system where [IC] is tightly con- 
trolled whereas [LDL] (the variable measured medically) is poorly controlled 
and can undergo wide variations. 



Based on experimental observations ((Goldstein and Brownl . I1977T ). we have 



modelled the regulation of the nuclear synthesis of LDL receptors as a simple 
reciprocal feedback. However, the qualitative features of the system are robust 
to our choice of functional form. The re gulation t erm i n Eq. (4) is a particular 
case of the standard inhibiting feedback ( Murray . 1990j ): c(l— 4>lr) / {q+[IC] n ), 



with n = 1 and q <C [IC] . Note, however, that the presence of q and n does not 
change the bounds (6)-(9). Therefore, the low and high cholesterol asymptotic 
regimes and the asymptotic response of the system to time-varying inputs 
(Section 4.2) remain unchanged. We have also checked numerically that the 
system still undergoes transitions between low and high cholesterol states; 
q and n have the effect of displacing the bifurcation points and modifying 
the extent of the coexistence region. In general, for large q the effect of the 
feedback is reduced, thus making the transition to the high cholesterol branch 
easier; for n > 1 the feedback is more efficient, thus enlarging the region of 
existence of the lower cholesterol branch. 
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4 Physiological Implications of the Model 



4-1 Parametric Effects of Medical and Genetic Conditions 



4-1.1 Tangier Disease and Reverse Cholesterol Transport 

Tangier disease is a deadly disease in which vanishin g reverse cholesterol 
trans port leads to accumulation of cholesterol in all cells (lYoung and Fielding , 
1999). In our model, the parameter die depends directly on reverse cholesterol 



transport. Note that only the bounds on [/C] eg are inversely proportional to 
die- This means that low values of die lead to diverging levels of intracellular 
cholesterol, while the plasma cholesterol levels saturate (Fig. 3). 



4-1.2 Lipase Enzyme Activity 

In a review article, Stein and Stein ( 2003f ) report that reduced enzyme activity 
can lead to reduced plasma cholesterol concentrations in mice and rabbits. 
Similar results might hold for humans. A conclusion of our sensitivity analysis 
(Table 3) is that variations in ki produce sharp changes in [LDL] without 
affecting the other variables significantly. Further numerics (not shown) show 
that decreasing ki decreases [LDL]. 



4-1.3 Familial Hypercholesterolemia 



Familial hypercholesterole mia (FH) is a genetic disease caused by mutat ions 



in the LDL receptor gene ( Goldstein and Brown . 1977 : Hobbs et al. . 1990j ). In 



the heterozygous case of FH, which affects about 1 in 500 people in generic 
populations and even higher ratios in more inclusive populations, cells express 
around half the normal number of functional receptors on their surface lead- 
ing to a two-fold increase in [LDL] . The homozygous case is much more rare 
(about 1 per million) but more severe. In this case, a negligible amount of 
functional LDL receptors lead to a dramatic rise of LDL plasma levels. These 
diseased individuals frequently dye of myocardial infarction or stroke in their 
teens. Table 4 present s a classification of FH gene mutations with their phys- 



iological implications ( Hobbs et al. . 1990h . We also show how the mutations 



can be represented in terms of the parameters of our model. 



Most of the medical observations in iHobbs et al. 
classes 3 and 4 



(1990) 



refer to mutation 
and homozygous. In terms of our model, het- 



both hetero- 

erozygous mutations 3 and 4 would halve di and d^. This produces roughly a 
two- fold (up to four-fold) increase in [LDL] , as deduced from the lower bound 
in (7) and Table 2. One of the homozygous mutations blocks the binding of 
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Mutation Class 


Physiological effect 


Affected Parameter 


l 


Reduced synthesis of LR 


ci 


Z 


Reduced transport of LR 


3 
4 


Reduced binding of IDL/LDL to LR 
Reduced number of LR in coated pits 


di i and/or di | 


5 


Diminished recycling of LR 


6T 



Table 4 

Gene mutations that le ad to famili al hypercholesterolemia. Mutations lead- 
ing to FH as classified bv lHobbs et alJ (|199(l ) and how they affect the parameters 
of the model. 

both IDL and LDL to LR, hence dj = di = 0. Again from (7), this would 
increase [LDL] by a factor (l + dL/d)(l + di/ki), which is roughly 10-fold (and 
up to several hundred t imes- fold) . These nu merical factors correspond well to 
those in the literature ( Hobbs et all Il990[ b Under another homozygous mu- 
tation 3 and 4, LDL receptors bind to IDL but not to LDL, i.e., di = 0. This 
produces an increase of [LDL] in the lower branch by a factor of (1 + di/d), 
which is two-fold (up to nine-fold) but still an order of magnitude smaller than 
that of the other homozygous case. Moreover, the characteristic time scale for 
the catabolism of LDL (see Section 4.2) increases from (rf + d^) -1 (around two 
days) to d^ 1 (over a week). Again, this coincides with experimental observa- 
tions that the removal of LDL from the plasm a is delayed and the increase in 
plasma LDL is moderate ([Hobbs et al. . 1990l b 



Medical evidence for mutations 1-2 and 5 is scarcely reported. In our model, 
they would translate into an increase of e = b/c which would have a small 
effect on the plasma concentrations on the lower cholesterol branch. However, 
as seen in Eqs. (15) and (16), an increase in e would lower the bifurcation 
point, making it easier to jump to the high cholesterol branch in response to 
an increase in uy or a decrease in die- Therefore, under these mutations the 
parameter region where the lower cholesterol regime exists would be reduced. 



Other observations in Hobbs et al. ( 1990j ) could also be interpreted in terms of 
our model. For example, it was reported that despite treatment with statins 
and lipid-lowering medications, the plasma cholesterol of some patients re- 
mained between 600 and 1200 mg/dl, five to ten times the normal levels. As 
seen in Figs. 2 and 3 and discussed further in Section 4.2.1, the presence of 
bistability means that once on the upper branch it is more difficult to return 
to the low cholesterol state by reducing %, e.g. through statin medication. 
Another intriguing observation is the dramatic clinical variability among FH 
patients with the same genetic mutation. A possible explanation is the ro- 
bustness of [/C] e q to variations in the physiological parameters. Indeed, from 
the lower bound in (9) we estimate that [IC] would only increase by 1.1 to 
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1.5 times under homozygotic FH, while the corresponding increase in [LDL] 
would range from one to almost three orders of magnitude. This would allow 
the development of individuals who would have similar levels of intracellular 
cholesterol while the LDL plasma levels could be very high in some cases, with 
different clinical consequences. 



4-2 Response to Time- dependent Inputs 



The time scales of the system span from a few hours up to a week. These 
time scales govern the way in which the system responds to external drives or 
perturbations to the variables. To clarify their dynamic al role, we use singular 
perturbation methods to obtain slow and fast manifolds (jKhahl koOO) relevant 
to the system. We have seen previously that the lower bound (4>lr = 1) and 
upper bound (4>lr = 0) are approached asymptotically as uy — > and as 
uy — > oo, respectively. In these singular asymptotic limits, the system (l)-(5) 
reduces to a linear system of the form: 



— = A k x + b, 
at 

where x = [[VLDL] [IDL] [LDL] [IC]] J 



{l,u}, 

= [Uy 0] T , 



(17) 
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-die 
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and l,u denote lower and upper asymptotic regimes, respectively. The inverse 
of the matrix eigenvalues are the characteristic times for each variable to return 
to equilibrium: 



(Jfej + diY 1 ss 0.5 < dj c w 1 h < K 1 



10 



< A;7 1 ^40< (d + ^ 70 < cT 1 ^ 200, (IS 



all of them given in hours. Note that the slowest time scale increases from 
(d + di)" 1 (on the scale of a day) on the lower branch to d^ 1 (over one week) 
on the upper branch. 



4-2.1 Response to slow drifts in the control parameters 

The response of the system to a slowly drifting input uy(t), which varies on 
a time scale slower than d" 1 , will be quasi-static, tracking the equilibrium 
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time (years) ' time (years) 

Fig. 4. Quasi-static nonlinear effects (a) Hysteresis curve obtained when uy is 
increased by 0.01 g{lh)~ l every 1000 hours from uy = lg{lh)~ l to uy = 2g(lh)~ 1 , 
and then decreased at the same rate to its initial value, (b) Bottleneck dynamics 
caused by the 'ghost' saddle-node around the bifurcation. The jump to the upper 
branch occurs a year after a small, but permanent, perturbation of uy has occurred. 



points. This can reveal nonl inear behaviour, such as hysteresis and bottleneck 
phenomena ( Strogat3 . 1994 b which can have physiological ramifications. 



If the system is in a parametric region where it can undergo bifurcations (as in 
Figs. 2 and 3), it will exhibit standard hysteresis when uy is swept continuously 
from the lower branch to the upper branch and backwards (Figure 4a). This 
means that once in the upper branch it is more difficult to return to the low 
cholesterol state by reducing uy (e.g., through statin medication). 

Another effect of the existence of saddle-node bifurcations is the possibility 
of bottleneck dynamics caused by saddle-node ghosts. This can translate into 
long-lived transients when the system is parametrically close to the bifurcation 
points. Figure 4b shows the time evolution of the system for two values of 
uy, one below and one above (but very close) to the bifurcation. The former 
remains in the low cholesterol branch and the latter jumps the high cholesterol 
solution. However, although the actual transition is relatively sudden (on the 
order of a week, ri _1 ), it occurs after a prolonged period of more than a year 
in which the system behaves as if in the (no longer existing) low cholesterol 
state. This pseudo-steady (ghost) behaviour would make it difficult to detect a 
very slow uy drift, which can suddenly lead to high plasma cholesterol values, 
if clinical measurements are infrequent. 
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4-2.2 Response to periodic inputs 




Fig. 5. Forced oscillations with a period of one day. Response of [IDL], 
[LDL] and [JC] to a periodic forcing uy = Uy (1 + sin(u;t)) with u = 27r/24/i~ 1 
for values of Uy = 1.25 (solid line), 1.75 (dashed line), 2.25 (dot-dashed line). All 
other parameters as in Fig. 2. The variables of the system exhibit different dynamical 
responses at this forcing frequency. The variability of [JC] is linked to that of [IDL] 
on the lower branch and to that of [LDL] on the upper branch. 



We study now the physiological implications of the response of the system to 
forced input oscillations with periods of days up to months. It is reasonable 
to assume that the production of VLDL (uy(t)) can have a strong periodic 
content on the order of a day since VLDL is produced at increased levels 
during the night; dietary intake is often highly cyclic; and the ingestion of 
statins to lower uy generally follows a daily cycle. As expected, entrainment 
phenomena can be observed in Figure 5, where we show the response of the 
system to forced oscillations with period of one day and increasing amplitude 
and mean. [LDL] is relatively insensitive to the oscillatory drive, while [IDL] 
exhibits large oscillatory amplitudes both for small and large drives. However, 
the salient feature is the different amplitude of the response of [IC] in both 
regimes. Qualitatively, this can be understood from the fact that [IC] changes 
from being tightly linked to [IDL] at low drives, when the amount of LDL 
receptors is high, to being governed by [LDL] on the upper branch, where the 
low (t>LR makes IDL endocytosis virtually impossible while LDL can still be 
absorbed through the non-receptor mediated route. 
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A quantitative understanding of the responses in Figure 5 emerges from a 
more detailed analysis. If we consider a periodic input 

u v (t) = U V (1 + Psm(cut)), (19) 

it follows that the time response of the asymptotic system (17) is sinusoidal 
of the form 

Xi = x i>av + x itS sm(ut) + x i>c cos(ujt), « = !,..., 4, (20) 



with average Xi >av and amplitude Xi tamp = ^xj s + xf c . We have calculated 

the sinusoidal responses for [IDL], [LDL] and [IC] on both branches. On the 
lower branch we get: 

[IDL Ul = -^-, = ■ P (21) 

k I + dl [wL Ul Jf^^yj J^TX? 



r rnrl kiUy [LDL] amP) i [IDL] amp j 1 
[LDL\ aVt i = — w r , r n = r n , == (22) 



[/C] „, 



(fcj + dj) (d + d L ) ' [LDL] aW) , [/DL] ot))J 

_ (x/dj + XLki) Uy 
(ki + dj) d/c 

[^C]amp,/ _ [IDL]amp,l 

[ICUj ~ [IDL] aVjl 



1 + 



\( 1 + (^k)' 

while on the upper branch we obtain: 



1 _|_ ( <£ Xi di 
~ \ d+d L Xldi+XLkr 



1 + te) 2 



(23) 



[JM]„„ = ^. ' rn r 1 "'' = , (24) 





[IDL] amp,u 




[IDL] aV)U 




[LDL] am p^ u 


d ' 


[LDL] aVjU 


XlU v 


[IC] amp,u 


die 


[IC] av ^ u 



2 

fe7 



L^ij — = — , rrnn = j^j— ^ (25) 

jj^fj _ /UA- y ^ J amp,u _ [LDL]amp,u 1 (26) 

die ' [IC] av ,u [LDL] av>u ^ ^ ^ 

The sinusoidal solutions (20)-(26), which are valid asymptotically in the limits 
uy — > (lower) and «y — > oo (upper), provide a good approximation when C/y 
is either small or large. Clearly, one does not expect a pure sinusoidal response 
at intermediate values of Uy, as can be observed in Figure 5. 

The change in the variability of [IC] on both branches can now be traced 
to Eqs. (21)-(26). For a forcing period of one day (uj = tt/12 0.26 h~ r ) 
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and using the time scales (18), the key terms are u/(d + dz) and uj/d. Note 
how on the lower branch (23), the u/(d + dz) terms cancel, with the effect 
that the relative variability [IC] amPt i/[IC] av> i is tied to that of [IDL]. Onn 
the upper branch (26), u /die ^ 1> which implies that the relative variability 
[lC\ amp , u /[IC]av,u is linked to that of [LDL]. 




time (years) d lc {h~ 1 ) 

Fig. 6. Forced oscillations with a period of 100 days, (a) Response of [LDL] 
to a periodic forcing uy = (1 + 0.2 sin(w t)) with u> = 27T/2400/1" 1 for values of 
die = [0.1, 0.2, 0.35, 0.45, 0.5], from top to bottom. All other parameters as in Fig. 3. 
As die increases, the system moves from the upper to the lower branch and the av- 
erage level of [LDL] decreases, (b) Relative amplitude of the oscillations, calculated 
as the ratio of range and mean of the oscillations from numerical simulations. The 
circles correspond to the values of djc in (a). Note that the relative amplitude is 
similar in both asymptotic regimes. This can be understood from Eqs. (22) and (25), 
whence we obtain the relative amplitude to be ps 0.38 (as die — * 0) and w 0.40 (as 
die ~~ * °°)- The relative amplitude reaches a maximum close to the bifurcation, 
where the system jumps between the upper and lower branches and is at its most 
nonlinear. As an indication, the bifurcation for these parameters and uy = 1 occurs 
at die ~ 0.373, close to the observed peak. 

These observations suggest that the correlation between intracellular and 
plasma cholesterol concentrations could be used as a marker to detect the 
high cholesterol state. This effect would be most pronounced in hepatic cells 
but might still be observable in other cells and detectable through a skin 
test. (This would depend on how much IDL can penetrate peripheral cells by 
leaking through the capillary endothelium into the interstitial fluid.) A less 
practical alternative would imply repeated time measurements of [IDL] to 
check for large amplitude oscillations. 
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We have also explored the response of the system to oscillations with lower 
frequency, as would be induced by seasonal cycles. Figure 6 shows the response 
of [LDL] to oscillatory inputs with period of 100 days at different values of 
die as shown in Fig. 3. Consistent with Eqs. (22) and (25), the average level 
of [LDL] on the upper branch is high, of order ~ gT 1 133, while on the 
lower branch it is much smaller, of order ~ kidj 1 (d + d^ 1 ~ 0.7. However, 
the relative amplitude of the oscillations of [LDL] in both asymptotic regimes 
is of the same order. It is precisely at the bifurcation point, when the system 
oscillates between the low and high cholesterol regimes, that the oscillations 
are larger in relative terms (Figure 6). This could be used to characterise the 
system if sufficiently frequent measurements are available. 



5 Discussion 

Understanding lipoprotein metabolism is important because high lipoprotein 
concentration in the plasma is a major risk factor for atherosclerosis, the most 
common cause of death in Western societies. We have introduced a dynamical 
model of lipoprotein metabolism derived by combining physiological cascading 
and cellular processes. The model goes beyond compartmental models in that 
we have included both the nonlinear absorption of IDL and LDL by cells, 
and a feedback mechanism by which the cell regulates the number of its LDL 
receptors based on the concentration of intracellular cholesterol IC. 

This low-dimensional, nonlinear model shows bistability and hysteresis be- 
tween a low and a high cholesterol state. Our bifurcation analysis of the key 
control parameters (uy and die) can be related to physiologically meaningful 
processes: diet and statin medication in the case of uy] reverse cholesterol 
transport through HDL in the case of die- We have checked that the observed 
behaviour is relatively robust with respect to all other parameters. An im- 
portant outcome of our sensitivity analysis is that the most robust feature in 
the low cholesterol state is the concentration of intracellular cholesterol [IC], 
while the plasma concentrations [IDL] and [LDL] can vary widely. This indi- 
cates that plasma cholesterol is not a tightly controlled variable in the system, 
implying that very different values of plasma [LDL] can be compatible with vi- 
able cellular function. It is conceivable, however, that the high levels of [LDL], 
innocuous at the cellular level, could lead to the initiation of disease-linked 
processes at the physiological level. 

We have also obtained estimates for the characteristic time scales governing 
the dynamics of the model in the low and high cholesterol states. These time 
scales account for the dynamical responses of the system to forced input os- 
cillations with very different periods. The model shows dynamical behaviour 
which could be used to diagnose the state of the system if repeated mea- 
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surements of plasma and/or cellular cholesterol with the correct timing and 
periodicity are performed, as indicated by Figs. 5 and 6. If measurements are 
too sparse, the onset of permanent high cholesterol levels might go undetected, 
with the possibility of a difficult return to lower levels due to hysteretic ef- 
fects. Measurements taken periodically on the order of hours or months could 
provide complementary clues to the state of the system. If we assume that 
periodic inputs to the system are important, measuring the relative variabil- 
ity of LDL and cellular cholesterol (possibly through skin tests) could lead to 
markers for the detection of the intrinsic regime of the system. 



Our dynamical model is directly derived from the underlying physiological 
and cellular processes. However, some of our approximations are crude and 
several directions for improvement could be pursued. The level of the plasma 
cholesterol on the high cholesterol branch ([VLDL] eq + [IDL] eq + [LDL] eq ) 
is unphysiologically high in the current model. This suggests the presence of 
other physiological mechanisms of elimination that hav e not been included in 
the m odel. Another extension would include time delays ((Dugard and Verriest , 
1998f ) to account for mass transfer (e.g., the delayed intake of IDL and LDL via 



LR- mediated endocytosis) and control mechanisms (e.g., the delayed genetic 
regulatory loop in the synthesis of LR). Furthermore, the model could be 
extended to include explicitly the regulatory and synthetic pathways involving 
the intracellular LDL receptors. It would also be important to distinguish 
between hepatic and non-hepatic cells, as hepatic cells play a more central 
role in lipoprotein metabolism. In particular, the rate of lipoprotein uptake 
in liver cells is different; increased uptake of LDL in liver cells lowers the 
production rate of VLDL (i.e., uy) through an additional regulating feedback; 
and statins reduce the inhibiting effect of [JC] on the synthesis of LR in liver 
cells only. 



Especially important would be to extend our description of reverse cholesterol 
transport to include [HDL] explicitly as a variable. Indeed, 20% of the total 
plasma cholesterol in healthy individuals (40mg/dl) is carried by HDL. Our 
preliminary work in this direction indicates that the number of nonnegative 
equilibria remains unchanged in the enlarged system that includes [HDL] , and 
the qualitative features of the system should hold. Because the relationship 
between plasma [LDL] and [HDL] is of major importance in clinical trials, 
this would be a crucial element to enhance the applicability of the model. 



Further modelling of lipoprotein metabolism could bring light to several open 
questions. To name but a few, it is still unknown why t he plasma cholesterol 

conce ntrations in adult humans are so unnecessarily high ((Goldstein and Brownl. 

1977 ): what the exact effects of statins in the metabolic system are ((Davidson and Jacobson . 



2001); and which of the processes could be best targeted to reduce plasma LDL 



concentrations. 
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Appendices 



A Derivation of the feedback term 



The feedback term in Eq. (4) follows from a more detailed (but still simpli- 
fied) derivation that includes the three compartments of intra-cellular LDL 
receptors (complexed, recycled, nucleus-synthesised) as shown in Figure A.l. 
The kinetic rate equations are: 



d[LR. 



compj 



dt 
d[LR r 



dt 

d[LR n 
dt 



(3a (dj[IDL] + d L [LDL\) LR - e[LR comp ] 
e(l - -f)[LR comp ] - 5[LR rec ] (1 - <p LR ) 



Cl 

[IC] 



- k[LR nuc ] - 5[LR n 



'I - 



'LR) 



(A.l) 
(A.2) 
(A.3) 



5 



dt 



-a (dj [IDL] + di [LDL]) <j) LR + -([LR rec ] + [LR nuc ]) (1 - <f) LR ) . 

(A.4) 



If we assume steady-state conditions: d[Li ^° mpl = d[L ^ recl = d[L ^" c] = 0, and 

^ r j-i in _. I a \ r_n t ] _ r _ / roi- 




] a(dl\DL]+d L [LDL]) 



Fig. A.l. Schematic description of the feedback term. An expanded view of 
the portion of Figure 1 that is related to the processes involved in the synthesis 
and recycling of LDL receptors inside the cell. If we assume steady state for the 
intracellular LR compartments, and that the turnover of nucleus-synthesised LR (k) 
is large compared to the probability of its insertion into the membrane (5), these 
processes lead to Equation (4). The factor f3 performs the appropriate dimensional 
conversions to account for the fact that 4>lr is a non-dimensional surface occupancy 
fraction. Note also that 7 is the fraction of non-recycled LDL receptors. All the other 
constants are rate constants with the appropriate units. 
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